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, Abstract. I use high resolution N-body/SPH simulations to model the new proper motion of 
' the Large Magellanic Cloud (LMC) within the Milky Way (MW) halo and investigate the effects 
of gravitational and hydrodynamical forces on the formation of the Magellanic Stream (MS). 
Both the LMC and the MW are fully self consistent galaxy models embedded in extended cuspy 
ACDM dark matter halos. 

I find that ram-pressure from a low density ionized halo is sufficient to remove a large amount 
. of gas from the LMC's disk forming a trailing Stream that extends more than 120 degrees from 
fS| ' the Cloud. Tidal forces elongate the satellite's disk but do not affect its vertical structure. No 
stars become unbound showing that tidal stripping is almost effectiveless. 

, Keywords, methods: n-body simulations, galaxies: Magellanic Clouds, galaxies: interactions, 
C^' galaxies: kinematics and dynamics 
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d . 1. Introduction 

Recent HST proper motion measurements of the Magellanic Clouds by Kallivayalil 



O 

O 



X 



et al. (20006a, b) and Piatek et al. (2007) indicate that they are presently moving at 
^ velocities substantially higher (almost 100 km s"'^) than those provided by previous 
observational studies (van der Marel et al. 2002, Kroupa & Bastian 1997). Such high 
, velocities {v = 378 km s~^ and v = 302 km s~^ for the LMC and SMC, respectively) 
■ are close to the escape velocity of the Milky Way and consistent with the hypothesis of 
, a first passage about the Galaxy (Besla et al. 2007). A single perigalactic passage has 
^-H serious implications for the origin of the Magellanic Stream. It definitely rules out the 
OO : tidal stripping hypothesis (Ruzicka, Theis & Palous 2008) since in this scenario the loss of 
' mass is primarily induced by tidal shocks suffered by satellites at the pericenters (Mayer 
^ , et al. 2006) and the Stream would not have time to form before the present time. In- 
deed kinematical data suggest that the Clouds are now just after a perigalactic passage. 
^ . On the other hand, ram- pressure scales as w^, where v is the relative velocity between 
satellites and the ambient medium. The high velocities of the Clouds could therefore 
compensate the effect of the reduced interaction time with the hot halo of the MW and 
hydrodynamical forces would play a determinant role in forming the Stream. 
In Mastropietro et al. (2005)1 (hereafter M05) we have performed high resolution N-body/SPH 



simulations to study the hydrodynamical and gravitational interaction between the LMC 



and the MW using orbital constraints by van der Marel et al. (2002) and a present time 



satellite velocity of 250 km s""'^. We found that, after two perigalactic passages, the com- 
bined effect of tidal forces and ram-pressure stripping can account for the majority of the 
LMC's internal features and for the formation of the MS. More in detail, ram-pressure 
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stripping of cold gas from the LMC's disk produces a Stream with morphology and kine- 
matics similar to the observed ones, while tidal stripping has longer time-scales and is 
not efficient in forming stellar debris, consistently with the lack of stars observed in the 
Stream. Nevertheless, at each pericentric passage the LMC suffers tidal heating which 
perturbs the overall structure of the satellite reducing the gravitational restoring force 
and therefore indirectly contributing to the loss of gas. 

The main objection to this model, in light of the new proper motion measurements of 
the LMC, is that the time spent by the LMC within the hot halo of the MW would be to 
short to cover the full extension of the Stream (more than 100 degree) by ram-pressure 
mechanisms (Besla et al. 2007). Moreover, hydrodynamical forces would affect a galaxy 
only weakly perturbed by the gravitational interaction and stripping would result more 
difficult. 

In this work I present the results of N-body/SPH simulations where the interaction 
between the MW and the LMC is modeled according to the new proper motion measure- 



ments of Kallivayalil et al. (2006a) 



2. Galaxy models 

The initial condition of the simulations are constructed using the technique described 
by Hernquist (1993) Both the MW and the LMC are multi-component galaxy models 
with a stellar and gaseous disk embedded in a spherical dark matter halo. The density 
profile of the NFW halo is adiabatically contracted due to baryonic cooling. Stars and 
cold gas in the disks follow the same exponential surface density profile. We also ex- 
plore the eventuality of an extended LMC gaseous disk. In this model the gaseous disk 
is characterized by an additional constant density layer which extends up to eight times 
the scale length of the exponential disk. The MW model comprises also a small stellar 
bulge and an extended low density (n = 2 x 10~^ cm~'^ within 150 kpc from the Galactic 
center and n = 8.5 x 10~^ cm~'^ at 50 kpc ) hot (T = 10^ K) halo in hydrostatic equilib- 
rium inside the Galactic potential (M05). The MW model, with virial mass 1O^^M0 and 
concentration c — 11, is similar to model Al of Klypin et al. (2002) while the structural 
parameters of the LMC are chosen in such a way that the resulting rotation curve re- 
semble that of a typical bulgeless late- type disk galaxy. In details, the satellite has virial 
mass 2.6 x IO^^A/q, concentration c = 9.5 and the same amount of mass in the stellar and 
gaseous disk component (~ IO^Mq). The Toomre's stability criterion is always satisfied 
and the parameter Q set equal to 1.5 and 2.0 at the disk scale radius in the different 
LMC models. 



3. Simulations 

I performed adiabatic simulations using GASOLINE, a parallel tree-code with multi- 
stepping (Wadsley et al. 2004). High resolution runs have 2.46 x 10^ particles, of which 
3.5 x 10^ are used for the disks and 5 x 10^ for the hot halo of the MW. The gravitational 
spline softening is set equal to 0.5 kpc for the dark and gaseous halos, and to 0.1 kpc for 
stars and gas in the disk and bulge components. 

In my best model the LMC approaches the MW on an unbound orbit (Fig. [T]) with 
initial Galactocentric distance of 400 kpc and velocity ^ 190 km s^^. After the peri- 
galactic passage (at ^ 40 kpc) the velocity decreases faster than for a ballistic orbit as 
a result of dynamical friction. The escape velocity at a given LMC position is indicated 
by a red curve in the right panel and calculated assuming a spherical unperturbed host 
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Figure 1. Orbit (left panel) and orbital velocity (right) of the LMC. Vertical lines represents 
the present time, just after the perigalacticon. Present values and values at the perigalacticon 
are indicated. 



potential. Due to the effects of dynamical friction, at late times the satellite lies on a 
nearly parabolic orbit. At the present time {t ^ 1.78 Gyr, vertical lines in the plots) it 
reaches a velocity of 378 km s"^ at 49 kpc from the Galactic center, in good agreement 
with the new proper motion measurements. 

In choosing the initial inclination of the LMC I made the approximation that it does not 
change during the interaction due to the effects precession or nutation of the disk plane. 
At the beginning of the simulation the disk moves almost face-on through the external 
medium, and ram-pressure affecting the whole disk perpendicularly. In proximity to the 
perigalactic passage the velocity vector changes rapidly and the angle between the satel- 
lite's disk and the proper motion is close to zero. At the present time the simulated disk 
has an inclination of about 30 degree with respect to the orbital motion (Kallivayalil 
et al. 2006a) and is indeed moving nearly edge-on through the external hot gas, with 
ram-pressure compressing its eastern side. 

Fig. [2] illustrates the present time distribution of stars and gas originating from the 
LMC's disk. The stellar disk becomes elongated while tidal debris start forming after 
the perigalactic passage. But all stars stay bound to the satellite. Tidal heating does not 
perturb significantly the vertical structure of the disk that remains thin and does not 
create a warp, unlike what observed in M05. Bar instability develops at the perigalacticon 
only in the case of Q ~ 1.5. 

Ram-pressure strips nearly 2 x IO^Mq of gas from the LMC's disk forming a continuous 
Stream that lies in a thin plane perpendicular to the disk of the MW and extending up 
to ~ 140 degrees from the LMC (Fig. [3]). The location of the Stream in the Southern 
Galactic hemisphere is comparable to the values of h and I provided by observations. 
Contrary to M05 there is no LMC's gas above the Galactic plane. In M05 the material 
lying in the Northern hemisphere is stripped from the satellite during the orbital period 
preceding the present one and the Stream forms a great polar circle around the Galaxy. 
The lack of gas at 6 > in Fig. [3] is not due to inefficient ram-pressure at early times, 
but to the fact that, in order to reproduce the current location and velocity of the LMC, 
the satellite enters the MW halo exactly at 6 = 0. 
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The morphology of the Stream does not change significantly adopting an extended 
gaseous disk model (top panel of Fig. [3]), except for the region at the head of the Stream, 
which appears broader (bottom). 

4. Conclusions 

I carried out high resolution gravitational/hydrodynamical simulations of the interac- 
tion between the LMC and the MW using the orbital parameters suggested by the new 
HST proper motion measurements. I find that ram-pressure stripping exerted by a ten- 
uous MW hot halo during a single perigalactic passage forms a Stream whose extension 
and location in the Sky are comparable to the observed ones. The stellar structure of the 
satellite is only marginally affected by tidal forces. 

The numerical simulations were performed on the 1 SGI-Altix 3700 Bx2 at the Uni- 
versity Observatory in Munich. This work was partly supported by the DFG Sonder- 
forschungsbereich 375 "Astro- Teilchenphysik" . 

References 

Besla, G., Kallivayalil, N., Hernquist, L., Robertson, B., Cox, T. J., van der Marel, R. P., & 

Alcock, C. 2007, ApJ, 668, 949 
Hernquist, L. 1993, ApJS, 86, 389 

Kallivayalil, N., van der Marel, R. P., Alcock, C, Axelrod, T., Cook, K. H., Drake, & A. J., 

Geha, M. 2006, ApJ, 638, 772 
Kallivayalil, N., van der Marel, R. P., & Alcock, C. 2006, ApJ, 652, 1213 
Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 
Kroupa, P., & Bastian, U. 1997, New Astronomy, 2, 77 

Mastropietro, C, Moore, B., Mayer, L., Wadsley, J., & Stadel, J. 2005, MNRAS, 363, 509 

Mayer, L., Mastropietro, C, Wadsley, J., Stadel, J., & Moore, B. 2006, MNRAS, 369, 1021 

Piatek, S., Pryor, C, & Olszewski, E. W. 2007, ArXiv e-prmts, a,rXiv:07 12A76A\ ArXiv e-prints, 

Ruzicka, A., Theis, C, & Palous, J. 2008, ArXiv e-pnnts, arXiv:0810.0968' 

van der Marel R. P., Alves D. R., Hardy E., & Suntzeff N. B. 2002, AJ, 124, 2639 

Wadsley, J. W., Stadel, J., & Quinn, T. 2004, New Astronomy, 9, 137 



High velocity LMC and formation of the MS 



123 




Galactic longitude 



Figure 2. Present time distribution of stars (top) and gas (bottom) from the LMC disk in 

Galactic coordinates. 
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Fieure 3. Polar Droiection of the simulated stream in Galactic coordinates. Both the Dure 



